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Abstract. We study the large deviation function for the entropy production rate in 
two driven one-dimensional systems: the asymmetric random walk on a discrete lattice 
and Brownian motion in a continuous periodic potential. We compare two approaches: 
the Donsker-Varadhan theory and the Freidlin-Wcntzcll theory. We show that the 
wings of the large deviation function are dominated by a single optimal trajectory: 
cither in forward (positive rate) or in backward direction (negative rate). The joining 
of both branches at zero entropy production implies a non-diffcrcntiability and thus the 
appearance of a "kink" . However, around zero entropy production many trajectories 
contribute and thus the kink is smeared out. 
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1. Introduction 



Driving a system away from thermal equilibrium implies a non-vanishing entropy 
production rate in the surrounding environment due to heat dissipation. For mesoscopic 
systems, the energy thus exchanged with the environment is a stochastic quantity, and 
so is the entropy s m produced along a single trajectory. This observation forms the basis 
of stochastic thermodynamics pQ. For a system driven into a non-equilibrium steady 
state, the long-time limit of the corresponding probability distribution p(s m , t) defines 
a large deviation function (LDF) 

h(a) = lim >P(^*) (1) 
v 1 t^oo (s m )t v ; 

with normalized entropy production rate 



a 



m . (2) 

Under very general assumptions, this large deviation function obeys the Gallavotti- 
Cohen symmetry [2] [3] 0] 

h(-a) = h{a) + a. (3) 

Experimental tests of this relation have been performed in a variety of systems, 
e.g., a single trapped colloidal particle [HI E] and a mechanical oscillator [7j. Other 
large deviation functions also play an important role for the study of non-equilibrium 
systems El QU [EE] . 
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Figure 1. The two systems studied: (a) the asymmetric random walk of a particle on 
a lattice with discrete positions. Different forward rate k + and backward rate k~ can 
be interpreted as an effective driving force / = ln(fc + /£;~). (b) The continuous analog 
is a particle moving in a periodic ring potential V(x) and driven by a force /. 



We study two paradigmatic models for driven systems sketched in Fig. [T] the 
asymmetric random walk on a discrete lattice and Brownian motion of a colloidal 
particle in a one-dimensional periodic potential. The latter system has been studied 
extensively in the context of the fluctuation-dissipation theorem EH El HS1 HE] • In 
order to obtain explicit expressions for Eq. Q, two main approaches are discussed in 
the literature (for a comprehensive review and references, see Ref. p2]). The first is the 
Donsker-Varadhan theory, stating that the scaled cumulant generating function 

a(X) = lim -ln(e Asm ) (4) 

t— >oo t 

is related to the LDF through the Fenchel-Legendre transformation 

h(a) = sup{Acr — a(X)/(s m )}. (5) 

A 

This method has been used previously to obtain the LDF for the two systems we study 
here: the asymmetric random walk [I] and the continuous periodic potential |18j . The 
second approach due to Freidlin and Wentzell [19J focusses on paths instead of the 
scaled cumulant generating function a(A). The probability of a path x(r) obeys a large 
deviation principle 

P[x(t)] ~ e - J[x(T)1/e (6) 

in the limit of small noise e — > 0. The large deviation function h(a) can then be obtained 
from the optimal trajectory x*(t) minimizing J[a;(r)] using the contraction principle. 
While the Donsker-Varadhan theory makes no assumptions on the noise strength and 
is thus more general, it is also less intuitive. The Freidlin- Wentzell theory, on the other 
hand, is more appealing to physical intuition since it allows a more direct insight. 

The purpose of this work is two-fold: First, we provide a detailed study comparing 
for two specific models the two approaches to the calculation of the large deviation 
function. Based on the Freidlin- Wentzell approach, the second purpose is to study 
the origin of a particular feature of the LDF: The emergence of a "kink" around zero 
entropy production, an abrupt albeit differentiable change of the slope [TS] . This 
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"kink" has been observed in distributions of the velocity (which is closely related to 
the entropy production) in models for molecular motors [20J, and experimentally for a 
tagged granular particle [21] . Of course, this "kink" does not indicate a phase transition. 
Dynamic phase transitions are indicated by a jump discontinuity of either the first or 
second derivative of the scaled cumulant generating function a(X). Such transitions have 
been found in interacting systems, e.g., for the totally asymmetric exclusion process [22J, 
the symmetric exclusion process [23], and for models with kinetic constraints |24j . 

Basically two theoretical explanations for the physical origin of the kink have been 
proposed so far. Pleimling and coworkers [251 126] have claimed that the kink is a 
generic feature of LDFs for the entropy production that follows quite generally from 
the Gallavotti-Cohen symmetry Eq. (|3]). But the relation Eq. ^ only restricts the 
anti-symmetric part, i.e., it allows us to write h(a) = h s (a) — a/2 with h s (—a) = h s (a). 
Clearly, from Eq. ^ no statement can be inferred about /i s , which is a non- universal 
function containing the details of the dynamics and the driving. On the other hand, the 
anti-symmetric part is universal and originates from breaking time-reversal symmetry. 
As an alternative explanation, Budini [27] has attributed the kink to "intermittency" , 

1. e., the random switching between different dynamic regimes. Here we show explicitly 
that the optimal trajectory minimizing the stochastic action indeed corresponds to two 
different regimes depending on the sign of a: forward (a > 0) or backward (a < 0). 
Considering only the optimal trajectories, therefore, implies a non-differentiability in 
the large deviation function at a = 0, which, however, is smeared out by fluctuations. 

2. Asymmetric random walk 

We first consider a particle on an infinite one-dimensional lattice. The particle can hop 
to the left with rate k~, or hop to the right with rate k + . A single trajectory of length t 
is characterized by the number of hops to the left, n~, and to the right, n + . Travelling 
the distance n = n + — n~, the entropy produced in the surrounding medium is 

k + 

s m (n + ,n") = n\n— = nf, (7) 

where 

(s m ) = f(k+-k-) = fk + (l-e~f) (8) 

is the mean entropy production rate. For k + = k~ the rates obey detailed balance and 
clearly no entropy is produced on average. 

2.1. Donsker-Varadhan theory 

In a seminal paper on the fluctuation theorem for stochastic dynamics [I], Lebowitz 
and Spohn have obtained the LDF for the asymmetric random walk (ARW) from the 
cumulant generating function. The normalized function h(a) can be cast in the form 

h{a) = (fz)' 1 |cosh(//2) + za[&rsinh(za) - f/2] - ^(za) 2 + lj (9) 
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Figure 2. Asymmetric random walk: (a) Large deviation function h(a) for two forces 
/ = 4 and / = 10 (solid lines). Also shown is the limiting form Eq. (12) (dashed 
lines), (b) The derivative becomes steeper at a 

forward (solid line) and backward (dashed line) fluxes [Eq. ( 18 )] vs 
zcr, where z = sinh(//2). 



for higher driving force /. (c) The 
the scaled rate 



with z = sinh(//2). It is straightforward to check that the function h(a) has the 
following properties: (i) it is convex and non-negative, h(a) ^ 0, with minimum 
h(l) = 0, (ii) it obeys the Gallavotti-Cohen symmetry Eq. (J3]), and (iii), following 
from this symmetry, h(—l) = 1 holds. 

To make our introductory comments concrete about a kink in the large deviation 
function, for \za\ ^> 1 we expand the inverse hyperbolic sine leading to 

h(a) « (fz)- 1 {cosh(//2) + za[±\nz\a\ - f /2] - z\a\} . (10) 

The plus sign holds for a > while the minus sign holds for o < 0. Here and in the 
following we assume / > 0. Considering the case / 3> 1 of large asymmetry we can 
further simplify 

C ° sh W 2 > » 1, (11) 
z ' 2 v ' 

yielding the limiting form 




l + a(lna-l) {a > 0) 

1 + -1 + f) (cr<0) 



(12) 



of the large deviation function in agreement with Ref. [25]. In particular, the first 
derivative 

A(ff)K 7i-InM-/ (*<0) (13) 

jumps at cr = 0. However, this jump is only apparent since for the expansion leading to 
the limiting expression Eq. (12) we used |er| 3> 1/z, i.e., we exclude a region of size 1/ z 
around a = 0. With increasing / this region becomes smaller but it does not vanish. 
Hence, for any finite / the function h(a) is differentiable. In Fig. [2j the large deviation 
function and its first derivative are plotted for different forces /. 
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2.2. Alternative derivation 

We now give an alternative derivation of Eq. ^ in the spirit of Freidlin-Wentzell. The 
path weight of a single trajectory is 

'> / " • ■■ • y~ e -(k++k~)t 

with m = n + + n~ total jumps. This weight fulfills the normalization condition 



P(n+,n-;t) = e ~^+* * (14) 



oo m 



^ £p(n+ «-;*) = 1. (15) 

m=0 n+=0 

Introducing dimensionless fluxes u ± = n ± /(k + t) and expanding the factorials using 
Stirling's approximation we obtain the rate function 

J(z/ + ,z/-) = - -InP 

= k + [l + e~ f - (u + + v~) + fv~ + u + \nu + + i/ _ lni/-], (16) 



which does not depend on time explicitly. For long times t the path weight Eq. (14) is 
sharply peaked around the minimum of J. We will call trajectories that contribute to 
this minimum "optimal" trajectories. Note that the small parameter e in Eq. Q is now 
given by e = 1/t, i.e., here the time limits the fluctuations around the minimum. 

For fixed a, forward and backward fluxes are not independent but related through 

u + = v- + (1 -e- f )a. (17) 

From d J / dv~ = we then find that the minimum of J is reached for fluxes 

\J{za) 2 + 1 ± za] , z = sinh(//2). (18) 



Using that 

t 

hxvf = — — ± arsinh(,2a~) (19) 

it is straightforward to show that J(u^,p~)/{s m ) = h(o) yields the result Eq. ^ 
obtained from the Donsker-Varadhan theory. 

The advantage of the alternative derivation is that it becomes transparent which 
trajectories contribute to h(a) for a given normalized entropy production rate a. We 
now use this insight to discuss the physical origin of the apparent kink. What does 
the optimal trajectory to reach a specified entropy production, i.e., the trajectory that 



minimizes the rate function J, look like? As becomes clear from Eq. (18), for a given 
positive a 3> 1/z the optimal trajectory consists almost entirely of forward jumps. 
Conversely, for negative a the optimal trajectory consists almost entirely of backward 
jumps even though these trajectories are extremely unlikely. This might seem surprising 
but any forward jump has to be compensated by an unlikely backward jump. Setting 
n~ = for a > and n + = for a < leads to the two branches of the limiting 



form given in Eq. (12). The reason why nevertheless h(a) is differentiable is that for 



\a\ < 1/z many trajectories contribute to the large deviation function. In particular, 
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for a = the forward and backward fluxes become equal, v + = v . All trajectories 
with the same number of forward and backward jumps contribute since they have the 



same weight. Their number is given by the binomial coefficient in Eq. (14). In contrast 



there is only one optimal combination (n + ,n ) in the wings of the distribution. 
3. Driven Brownian motion 

In this section we consider a single colloidal particle moving in a closed ring geometry 
with circumference L and external periodic potential V(x+L) = V(x). The overdamped 
dynamics are described by the Langevin equation 

x = fi [-V'(x) + f}+ V , (20) 

where the prime denotes the derivative with respect to x. The particle is driven out of 
equilibrium through the force / resulting in a non-vanishing current through the ring. 
The noise r] describes the effective interactions between the solvated particle and the 
solvent at temperature T with correlations 

(rj(r) V (r')) = 2D 5(r - r'). (21) 

Through Dq = [i^T we assume that the solvent remains at equilibrium even though the 
particle is driven. To simplify the notation we will use dimensionless quantities 

x^Lx, V(x)^TV(x), f^yf, r^^-r (22) 

throughout the remainder of this paper. Moreover, these units make transparent the 
connection with the ARW. In particular, the number of periods traversed is n and the 
dimensionless force / corresponds to the force defined in Eq. ([7]). We will also need the 
notion of a critical force f c such that for f > f c the deterministic force —V'(x) + / > 
is positive for all x and, therefore, in the absence of noise running solutions exist. 

The rate of entropy production is s m = fx with s m = f Ax [cf. Eq. (m) for the 
ARW], where Ax is the distance the particle has traveled during the time t. The mean 
velocity of the particle can be calculated from Stratonovich's formula [28] 

(£) = k(1 - e -/ ) (23) 

with 

k ^ = [ dx [ dz exp{V(x) - V(x-z) - fz) . (24) 
Jo Jo 

The mean velocity (and also the mean entropy production) of ARW and ring potential 
agree if we set k + = k, cf. Eq. ([8]). However, in general fluctuations will be different, 
leading to different shapes for the LDFs. For large forces / the influence of the potential 
on the particle motion becomes negligible and (x) ~ /. 
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3.1. Donsker-Varadhan theory 



The time evolution of the probability distribution of x is governed through the Fokker- 
Planck operator 



L n = 



d_ 
dx 



(V -f) + 



d_ 
dx 



In particular, the steady state distribution obeys Lop ss (x) 
[Eq. Q] consider the two-sided Laplace transform 

" + OC 



p(x,X,t) 



ds r 



l p(x,s m ,t) 



(25) 

0. To calculate a (A) 
(26) 



of the joint probability p(x,s m ,t) to find the particle at position x with an amount 
of entropy s m produced up to time t. The time evolution of the transformed joint 
probability obeys d t p = L\p, where the operator 

d 



L + 



2 ir + v ' 

ox 



f 



(/a) + (f\y 



(27) 



has been determined in Ref. 
the time evolution operator 



. We expand the joint probability in eigenf unctions of 



P(x,\,t) = j2^ ai{x)t ^( x > x )' 



(28) 



8=0 



where the coefficients c, quantify the overlap with the initial distribution p(x, X, 0) = 
p BB (x). The eigenvalue equations read L\ipi(x,X) = ai(X)ipi(x, A). In the limit of long 
times, the joint probability is dominated by the largest eigenvalue a . Integrating 
Eq. (26) over the position x thus leads to 



and we conclude that the scaled cumulant generating function a (A) 



(29) 

ao(A) is indeed 



given by the largest eigenvalue of the time evolution operator (27). It can be calculated 
solving either the eigenvalue problem directly [18J or through the Ritz variational 
method [29]. The function h(a) then follows from Here we use the first method 
and refer to Ref. HHI for further details. 



3.2. Stochastic action and effective potential 

To continue we employ the well-established machinery of path integrals for stochastic 
processes (30] that has been developed to deal, e.g., with diffusion in bistable 
potentials [31] and escape problems [32] . Imagine that the particle starts at position x$ 
and travels a distance Ax during the time t. The probability for such a transition 

i>x +Ax 

P{x + Ax\x ;t) = / [dx{T)}e- s[x{T)] (30) 

•Ixo 

is obtained through summing over all possible paths with stochastic action 

S[x(t)]=J\t ^{x + V'-fY- l -V"} (31) 
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calculated along a single path x(t). The last term V"/2 is the Jacobian arising from 
the change of variables 77 — > x. 

The stochastic action Eq. (31) can be rearranged to give 



1 . 1 „ . 1 

6>(r)j = 

The last term 



S[x{r)\ = -AV - -fAx + -ft + I dr C(t). (32) 







C(x,x) = he 2 - U(x) (33) 

takes on the form of a Lagrangian for a classical particle with mass m— 1/2 moving in 
the effective potential 

U(x) = l -V"{x) - \[V\x)f + \fV\x). (34) 

It is straightforward to check that for a periodic (and twice differentiable) potential V(x) 
the effective potential U(x) is periodic, too. For completeness, we note that the same 
effective potential arises in the symmetrization of the Fokker-Planck operator [Eq. (25)], 

e y^L e- v ^ 2 = ^- 2 +U(x), (35) 



see, e.g., the textbook treatment in Ref. 

In Fig. [3] the effective potential U(x) is plotted for the cosine potential V(x) = 
t> cos(27rx) with critical force f c = 2tcv . It will turn out useful to define the reduced 
force /* = f / f c - In the following we will use this specific potential for all illustrations. 
In equilibrium, the principal feature of the effective potential is that the minimum of 
V(x) becomes the global maximum of U(x). For v > 1, the minimum of U(x) splits into 
two minima, which are degenerate at equilibrium (/ = 0). Turning on the driving force 
(/ > 0) the global maximum of U(x) is shifted, the depth of the left minimum increases, 
and the depth of the right minimum decreases until for forces /* ^ (1 — v 2 ^ 3 ) 3 / 2 only 
one global minimum remains. For /* ^ 1 the upper arc becomes concave. 




Figure 3. (a) Original potential V(x) = vq cos(27ra;) and (b) the mapped effective 
potential U(x) from Eq. (34) plotted for vq — 3 and forces / = 0, / < / c , and at the 
critical force / = f c = 6ir. 
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3.3. Optimal trajectory 



We now determine the trajectory x*(r) that minimizes the action [Eq. (32)] with 
<S* = <S[x*(t)]. Clearly, this trajectory is the solution of the Euler-Lagrange equation 

5S _ d dC dC 1 „^ u'(x) — (36) 
8x(t) dr dx dx 2 

subject to the boundary conditions x*(0) = x and x*(t) = x + Ax. Since the effective 
potential is periodic, the solution x*(t + t ) = x*(t) of the Euler-Lagrange Eq. (36) will 
also be periodic with period 

1 



To 



\*\{x) 



(37) 



The sign of a determines the sign of the initial velocity of the particle. Note that the 
solution of Eq. (36) is time-reversal symmetric. 



Before we proceed we shift the effective potential 
U(x) = U* + W(x), U* = max U(x) 



(3* 



such that W(x) ^ is non-positive. Although no small parameter appears in the 
stochastic action, we nevertheless follow the Freidlin-Wentzell route and assume that 
the transition probability is dominated by the optimal trajectory. We thus obtain from 
Eq. (32) the unnormalized, naive approximation 

<5* 



h. 



lim 

t— ¥O0 t 



1 /• V 1 s-2 A 



u. 



with optimal stochastic action for a single period 



A = 



dr 



[x*(t)] 2 -W(x*(t)) 



(39) 



(40) 



Since the potential V(x) is bounded, the term AV/t drops out in the long time limit. 
We stress that, since no small parameter appears in the stochastic action, it cannot be 
expected in general that h*(a) corresponds to h(a) as obtained by the operator method. 

The time integral in Eq. (40) is evaluated most conveniently by changing the 
integration variable from r to x. To this end we note that 
1, 



E 



is a constant of motion. The time to traverse a single period is 

2 Jo 



ro(E) 



dx [E - W(x)\ 



-1/2 



Rearranging Eq. (40) using Eq. (41), the action for a single period reads 

-I 



A(E) 



dx 



o 



E - 2W(x) 
2JE - W(x) 



dx y/E-W(x)-ET (E). 



(41) 



(42) 



(43) 



o 



We thus solve Eq. (|42|) for a given time To in order to obtain E, from which we determine 

(44) 



the optimal action A. The resulting normalized function 
K{a) = K(a)/(s m ) 




Figure 4. Large deviation function for the periodic potential V(x) = vo cos(27rx) with 
vq = 2 (top row) and vq = 6 (bottom row) for three reduced driving forces /* . Shown 
are the exact solutions h(cr) (thick lines) obtained from the operator method together 
with the naive approximations /i*(cr) [Eq. (44), thin lines]. For small a the numerical 
calculation breaks down and the dashed lines show the linear extrapolation towards 
the known value h*(Q) [see Eq. (49)]. 



is plotted in Fig. [4] for different driving forces / below and above the critical force f c . 
The calculation breaks down for small a as the initial velocity, and therefore E, become 
of the order of the machine precision. We find that for large a the optimal path solution 
agrees with the full solution h(a) but deviates for small a, where the disagreement 
becomes more pronounced for small driving forces. 

We now consider two special cases. For large E the influence of the potential is 
negligible and the velocity ~ ±2yE is approximately constant, with x* = ±l/r = 



a(x). Together with A = Etq we obtain from Eq. (39) 

K(<T) = ±(<r(x)-f) 2 -U<, (45) 

i.e., the asymptotic form of the large deviation function for large rates |cr| is a parabola. 
In particular, for a flat potential V(x) = with (x) = f we recover the exact result 
h{a) = {l/A){a-lf. 

For E — > the period Tq increases dramatically and the particle spends a long time 
at the barriers of the effective potential [corresponding to the valleys of the original 



potential V(x)], see Fig. [5j Up to linear order in E the optimal action Eq. (43) becomes 

A ~ A = A(0) = f dx y/-W(x). (46) 
Jo 

This is of course the action of an instanton, the classical solution for E = that 
connects two adjacent barriers in infinite time. To approximately calculate the period 
T we expand W(x) to second order, 



2 ~ 2 



dx [E - [7^72] - 1 / 2 , (47) 
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Figure 5. Optimal trajectory x^(t) plotted over three periods (parameters: Vq = 5, 
/ = 8tt, a = 1). The potentials U(x) (red/bright) and V(x) — fx (black) are shown in 
the right panel. 



where, without loss of generality, we have shifted the potential along the x-axis such 
that W(0) = 0, and U£ < is the curvature of the effective potential at the barrier. To 
leading order we find the well-known result [31] 

£^M exp j_y^ ro j. (48) 

While the time r to traverse a single period diverges for a — > 0, the action of an 
instanton remains finite and Aq/tq — > 0. The approximated large deviation function at 
vanishing entropy production a = thus reads 

~K(U) = \f 2 -U*. (49) 

The linear extrapolation of the numerically determined h*(a) towards this value is shown 
in Fig. |4| 

Summarizing, the results of naively applying the Freidlin-Wentzell approach to the 
calculation of the large deviation function are: (i) The optimal path is not sufficient 
to reproduce the large deviation function over the full range of rates a. (ii) Only 
accounting for the contribution from the optimal path results in a function h*(a) that 
is non-differentiable at a = 0, i.e., the approximated large deviation function exhibits 
a genuine kink, (iii) The Gallavotti-Cohen symmetry ^ is fulfilled by both h*(a) and 
h(a). The rounding of the kink in the full solution h(a) points to the importance of 
fluctuations around the optimal path even in the limit t — )■ oo. 



3.4- Small-noise limit 

As mentioned in the introduction, the Freidlin-Wentzell method is strictly valid only 
in the limit e — > of vanishing noise. In our dimensionless units, the small parameter 
e = 1/vo is the ratio of thermal energy to the potential depth. Fixing /* with f c oc vq 
we see that the first term in Eq. ( 34 ) is of order vq and the other two terms are of order 
v^. We thus drop the first term (the Jacobian). Rearranging terms we obtain for the 
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Figure 6. Small-noise limit: (a) Simplified effective potential U(x) dropping the 
Jacobian in Eq. ([34} [cf. Fig. ^b)]. (b) Approximated LDF h*(a) for /* < 1 (thick 
line). Also shown is the LDF h(a) for vq — 7 and / = 6ir (/* = 3/7) together with the 
ARW solution Eq. ^ for the same force. 



effective potential the simplified expression 

U(x) = \f - \[-V'{x) + ff. (50) 

In the following we focus on the case / < f c . The maximum of the potential is then 
reached at points x where the deterministic force —V'(x) + / is zero with 

U. = \f (/ < f a ). (51) 



From Eq. (49) we then know that ft*(0) = 0. Since h*(a) is a convex function with 
/i*(l) = we can infer that the shape in the range — 1 ^ a ^ 1 is a "wedge". The 
approximated LDF is universal in this range and does not depend on / (as long as 
/ < f c ) nor vq. 

Specifically, for the cosine potential V(x) = vq cos(27r:r) the effective potential 

U(x) = & {^-[sin(27rx) + /,] 2 } (52) 



is shown in Fig. [6[a), whereas the "wedge" -shaped h*(a) is depicted in Fig. [6[b). The 
action of an instanton is easily evaluated to give 

A = [ dx |sin(27rx) + /„| = — a/ 1 - / 2 + /* arcsin/* . (53) 
J 2 it L J 



Also shown in Fig. |6[b) is that the actual LDF h(a) agrees with the ARW solution for 
large vq. Since /* is fixed the force / increases as we increase the potential depth u . 
Both the ARW and the actual solution h(a) then approach the limiting "wedge" shape 

as Vq —> OO. 



3.5. Role of fluctuations 

At this point it seems natural to inquire what the nature is of the fluctuations 
contributing to the large deviation function. Specifically, we consider two types of 
fluctuations: small Gaussian perturbations around the optimal path, and fluctuations 
of the "jumps" times 77. In the periodic optimal trajectory x*(r) these jumps occur 
regularly with period r , see Fig. [5j 
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3.5.1. Gaussian fluctuations. For finite t, we expand the action Eq. (32) 

S[x.{r) + £(t)) ^S*+ l -j\r i{r)D{r)i{r) (54) 

to second order with small perturbations £(r) of the optimal path, which obey the 
boundary conditions £(0) = = 0. Here we have introduced the symmetric second 
variation operator 

£>(r) = -~-U"(x.(r)), (55) 

which is formally equivalent to a Schrodinger equation with potential energy —U"{x*(t)) 
(positive curvatures of the effective potential correspond to low "energies"). We, 
therefore, can immediately state that the eigenvalues Aj of D are real and can be ordered, 
Ao < Ai < • • - . Moreover, eigenfunctions £i(r) corresponding to different eigenvalues \ 
are orthogonal, 



Jo 



6(0) = &(t) = 0. (56) 



We expand the generic perturbation £(r) in the basis spanned by the normalized 
eigenfunctions {£.;}, 

oo 

£(r) = 5>&(r), (57) 
i=i 

with coefficients q. The path measure is 

oo 

[dx(r)] oc Y[dc t (58) 



i=l 



and we perform the Gaussian integration of Eq. (30) with result 

e~ s * 

P(x + Ax\x ; t) ~ — = , detD = j|Aj. (59) 

V det t) i= i 

Of course, this result only holds if all eigenvalues are positive. We numerically calculate 
the determinant using Floquet theory as outlined in the appendix. We find that small 
Gaussian fluctuation around the optimal path x*(r) do not contribute to the large 
deviation function. 



3.5.2. Dilute instanton gas. For / <C / c the particle in the periodic continuous 
potential V(x) effectively undergoes a hopping motion resembling the discrete random 
walk. For such thermally activated jumps, the corresponding rate is 

k + ^exp{-AV + f(x + -x-)} , (60) 

where x~ is the position of the minimum, x + the position of the barrier, and AV = 
V(x + ) — V(x~) is the barrier height. For the cosine potential we find 

2irx~ = 7T + arcsin /*, 2irx + = 2n — arcsin /*. (61) 
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The rate thus reads 
k + «< 



-A0+//2 



(62) 



where A is the action [Eq. (53)] of an instanton in the small noise limit and we have 
ignored a sub-exponential prefactor. 

The optimal trajectory x*(r) is strictly periodic. However, Fig. [5] suggests that 
shifting the positions of the "jumps" n = Itq does only minimally change the action as 
long as the particle spends a quasi-infinite time on top of the barrier (of the effective 
potential U). In this case some eigenvalues \ become exponentially small (and therefore 
contribute to the large deviation function) or even vanish, in which case the calculation 
of the previous subsection breaks down. Such zero-modes need to be treated in a 
different way known as the method of collective coordinates. In effect, integrations over 
zero-modes are replaced by integrations over the corresponding jump time. 

Suppose a trajectory is composed of n + instantons and n~ anti-instantons jumping 
against the driving force. We use again n = n + — n~ and m = n + + rT . If we neglect 
all "interactions" between instantons we can approximate the action as 

Szi-^fn + mAo (63) 

independent of the actual jump times. Integration over the jump times then leads to a 
path weight 



1 

Z 



m 

n + 



(64) 



P(n + ,n-;t) 

with normalization factor 

Z = exp {2t cosh(f/2)e- Ao } . (65) 

We now follow the alternative derivation for the asymmetric random walk presented in 
Sec. 2.2 Introducing dimensionless fluxes v ± = n^/fat) [related by Eq. (17)], we obtain 
the rate function 

J(v + ,v- 



2cosh(//2)e- Ao + k{v a 
_(//2)(i/+-0-(^ 



In v + + v In v 
+ u~)(l -A - Iuk)} 



The minimum is attained at fluxes 

,-A 



with 



v. 



c 



± 



Is' 



V(C^) 2 + 1 ± 



-ne 



A01 



(66) 



(67) 



(68) 



Plugging in the rate Eq. (62) for k = k + we find ( ~ z and, moreover, we recover the 
form Eq. (191) of the LDF for the ARW. 
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4. Conclusions 

We have studied the large deviation function for the entropy production in two simple 
one-dimensional systems: the asymmetric random walk of a particle on a discrete lattice 
and the continuous motion of a driven particle in an external potential. For both systems 
we have calculated the large deviation function using the path approach of Freidlin 
and Wentzell. However, while for the ARW the solution thus obtained agrees with the 
Donsker-Varadhan theory, for the continuous case the solution is only an approximation. 
In both approaches, the wings of the large deviation function are well described by a 
single trajectory composed of only forward (cr > 0) or only backward (a < 0) jumps. 
Changing the prescribed rate a, the transition between these two regimes appears as a 
kink at o = 0, which in the true large deviation function is smeared out by fluctuations. 
For the ARW the reason is that there are many combinations of forward and backward 
jumps (and therefore many trajectories) around a = as described by the combinatorial 



factor in Eq. (14). For the continuous system we have verified numerically that small 
perturbations of the optimal path do not contribute to the large deviation function. For 
large vq and small driving forces /, < 1 trajectories are dominated by discrete barrier 
crossings and can be decomposed into instanton solutions. Summing these trajectories, 
we recover the large deviation function of the ARW. 
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Appendix A. Floquet theory 

We estimate the effect of small fluctuations through calculating the full determinant 



det D of the operator Eq. (55). To this end we make use of the following relation 



detD _ <j>{t) 
det Do Mt) 

with D a suitable reference operator [33] • In the simplest case D = —d\ but the exact 
form of Do does not play a role. Here, <f>(t) is the solution of Dcf) = whereas (f>o(t) is 
the solution of D <f)Q = 0. Defining the vector x = (<f),<f)) we obtain a linear, periodic 
differential equation 

x = M(^(t))x, M(x) = ^ _J„ {x) I j . (A.2) 

From Floquet theory it is well-known that the solution can be written 

x(r) = c + e +fMT p{r) + c_e^ T p(-r) (A.3) 



CONTENTS 



17 



with characteristic, or Floquet, exponent fj,. Hence, in the limit of long times the 
determinant behaves as det D ~ <p(t) ~ |e M *| to leading order. We again change variables 
from r to x and solve the ordinary differential equation 

^ = —M(x)X = \[E - U(x)}- 1/2 M(x)X 



for a matrix X with the matrix M(x) given in Eq. (A.2). The initial condition is X(0) = 1, 
and we integrate a single period. The resulting matrix X(l) has two eigenvalues 
p± = e ±At . Numerically we find \p±\ ~ 1 for the values of E that can be obtained 
numerically. Therefore, small Gaussian fluctuations around the optimal path x*(r) do 
not contribute to the large deviation function in the long time limit. 
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